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The dynamics of infinite, asymptotically uniform, distributions of purely self-gravitating particles 
in one spatial dimension provides a simple and interesting toy model for the analogous three di- 
mensional problem treated in cosmology. In this article we focus on a limitation of such models as 
they have been treated so far in the literature: the force, as it has been specified, is well defined in 
infinite point distributions only if there is a centre of symmetry (i.e. the definition requires explicitly 
the breaking of statistical translational invariance). The problem arises because naive background 
subtraction (due to expansion, or by "Jeans swindle" for the static case), applied as in three dimen- 
sions, leaves an unregulated contribution to the force due to surface mass fluctuations. Following a 
discussion by Kiessling of the Jeans swindle in three dimensions, we show that the problem may be 
resolved by defining the force in infinite point distributions as the limit of an exponentially screened 
pair interaction. We show explicitly that this prescription gives a well defined (finite) force act- 
ing on particles in a class of perturbed infinite lattices, which are the point processes relevant to 
cosmological A'^-body simulations. For identical particles the dynamics of the simplest toy model 
(without expansion) is equivalent to that of an infinite set of points with inverted harmonic oscil- 
lator potentials which bounce elastically when they collide. We discuss and compare with previous 
results in the literature, and present new results for the specific case of this simplest (static) model 
starting from "shuffled lattice" initial conditions. These show qualitative properties of the evolution 
(notably its "self-similarity") like those in the analogous simulations in three dimensions, which in 
turn resemble those in the expanding universe. 

PACS numbers; 



I. INTRODUCTION 

The development of clustering in initiaUy quasi- 
uniform infinite distributions of point particles evolving 
purely under their Newtonian self-gravity has been the 
subject of extensive numerical study in cosmology over 
the last several decades (see e.g. jlj for a review). This 
is the case because these "N-body" (particle) simulations 
of the Newtonian limit are believed to give a very good 
approximation to the formation of structure formation in 
current dark matter dominated models of the universe. 
The impressive growth in the size of these simulations has 
led essentially to phenomenological models of the associ- 
ated dynamics. Analytical understanding, which would 
be very useful in trying to extend the numerical results 
and also control for their reliability, remains very limited. 
In attempts to progress in this direction it is natural to 
look to simplified toy models which may provide insight 
and qualitative understanding. Such models may also be 
interesting theoretically in a purely statistical mechan- 
ics setting, and specifically in the context of the inves- 
tigation of out of equilibrium dynamics of systems with 
long-range interactions (see e.g. 

An obvious toy model for this full 3 — d problem is the 
analogous problem in 1 — d, i.e., the generalization to an 
infinite space (static or expanding) of the so-called "sheet 
model" , which is formulated for finite mass distributions. 
In this latter model, which has been quite extensively 



investigated (see, e.g., |^, ^, |[ 0, ||, D and references 
therein), particles in 1 — d experience pair forces inde- 
pendent of their separation, like those between parallel 
self- gravitating sheets in 3 — d of infinite extent. Several 
groups of authors |lO[ |ri|, g Q ^, ^ 

have then discussed different variants on this model to 
develop the analogy with the 3 — d infinite space prob- 
lem. Just as for the finite sheet model, these models 
have the particular interest of admitting exact solutions 
between sheet crossings, which means that they can be 
easily evolved numerically to machine precision, and at 
modest numerical cost for quite large numbers of parti- 
cles. 

In this article we revisit the basics of these toy models 
(in either static or expanding universes), addressing the 
problem of their general formulation for infinite distri- 
butions. Indeed, as we will discuss, previous discussions 
have required, in their implementation, the imposition of 
symmetry about a point, or finite extent of the consid- 
ered density perturbations^. Such a restriction on the 
class of point processes which can be considered, and no- 
tably the requirement that statistical translational invari- 



^ This is not true of the treatments in , which start directly 

from the fluid limit (rather than from a particle description) . See 
further discussion below. 
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ance be broken, is not desirable. Indeed in the context 
of the cosmological problem, this latter property of the 
distributions usually considered as initial conditions for 
simulations is very important, because of the "cosmologi- 
cal principle" which supposes that there are no preferred 
centres (see e.g. §). Further the question of the 
extrapolation of the finite version of the model (which 
is what is simulated numerically) to the infinite system 
limit has, as we will discuss below, not been carefully ex- 
amined. We will show that problems with the definition 
of the force (as used in these previous treatments) arise 
from a subtlety about how the so-called "Jeans swindle" 
is applied in one dimension. We draw here on the work of 
Kiessling in |2^, where it has been shown that, in 3 — d, 
the usual formulation of the Jeans swindle — subtraction 
of a compensating negative mass background in calcula- 
tion of the potential — may be more physically formu- 
lated as a prescription for the calculation of the force in 
the infinite volume limit. It turns out, as we will see, that 
while in 3 — d it is sufficient to prescribe that the force on 
a given particle is obtained by summing symmetrically 
about it (e.g. summing in spheres of radius R with cen- 
tre at the particle, and then sending R to infinity), in 
1 — d this limiting procedure needs to be further speci- 
fied. More specifically the force turns out to be defined 
in 1 — d for a broader class of point distributions — and 
notably for distributions without a centre — when the 
summation is performed by taking the unscreened limit 
of the same sum for a screened version of the interaction, 
rather than as the limit of the sum truncated to a finite 
symmetric "top-hat" interval. 

In the next section we give a more detailed heuris- 
tic discussion of the problem of defining the force in the 
infinite volume limit, and then giv e the prescription we 
adopt. We then present in Sect. HI a rigorous calculation 
showing that the force is indeed well defined for a certain 
class of infinite perturbed lattices, i.e., infinite configu- 
rations generated by perturbing particles off a perfect 
lattice. To do this we treat these infinite point distribu- 
tions as stochastic point processes and study the prob- 
ability distribution of the force on particles p^, p4[ . For 
finite variance displacements off the lattice, either cor- 
related or uncorrelated, which do not cause particles to 
cross one another, the result turns out to be extremely 
simple: the force on any particle is simply proportional 
to its displacement from its lattice site. In the following 
section we turn to consider the definition of dynamical 
models using this force. While the most evident model is 
the simple one obtained by using the conservative New- 
tonian dynamics under the derived forces, there is also a 
simple variant with an additional damping term which is 
the natural toy model for the 3 — d cosmological prob- 
lem (with an expanding background). Given that par- 
ticle crossings are, up to interchanges of particle labels, 
equivalent, in 1 — d, to an elastic collision (with exchange 
of velocities) the evolution in these toy models, starting 
from infinite perturbed lattices in the class in which we 
have shown the force to be defined, is in fact equivalent 



at all times simply to an infinite set of particles with 
inverted harmonic oscillator potentials centred at their 
original lattice sites, and which collide elastically when 
they meet. In the last subsection of this section, we then 
give a detailed discussion of the relation of these two toy 
models to those which have been discussed previously in 
the literature, explaining that our formulation provides 
essentially both a simplification, and a generalization, of 
most of these previous treatments. In Sect. ^ we present 
briefly results of numerical simulations for the simplest 
toy model (without expansion), and uncorrelated initial 
displacements to the lattice (the "shuffied lattice"). We 
show that in this case the evolution of the clustering 
in time is qualitatively very similar to that which has 
been observed in the analogous 3 — d system. Notably 
these static space simulations share the features of "hier- 
achichal" structure formation and "self-similarity" which 
are well documented in full 3~d simulations. In the final 
section we summarize our findings and conclusions, and 
discuss some directions we envisage for further work. 



II. FROM FINITE TO INFINITE SYSTEMS 

A. Definitions 

By gravity in one dimension we mean the pair interac- 
tion corresponding to an attractive force independent of 
separation, i.e., the force f{x) on a particle at coordinate 
position X exerted by a particle at the origin is given by 



/(^) 



-g sgn{x) 



(1) 



where g is the coupling. Equivalently it is the pair in- 
teraction given by the pair potential (j){x) — g\x\ which 
satisfies the 1 — d Poisson equation for a point source, 

= 2gSD{x) (where Sjj is the Dirac delta function). 
Comparing with the 3 — d Poisson equation shows the 
equivalence with the case of an infinitely thin plane of in- 
finite extent and surface mass density S ~ g/2TrG, which 
explains the widely used name "sheet model". We will 
work in the one dimensional language, referring to "par- 
ticles". For convenience we will set the mass of these 
particles, which will always be equal here, to unity. 



B. Finite system 

Let us consider first the case of a finite system, con- 
sisting of a finite number N of particles (with either open 
boundary conditions, or contained in a finite box). De- 
noting by Xi the coordinate position of the i*'' particle 
along the real axis, the force field F{x) (i.e. the force on 
a test particle) at the point x is 

n.)^»y:.gn(..-x,^J..„(„.,gn(.-.),,2) 
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where n(y) = J^i ^oiv — Xi) is the microscopic number 
density and the integral is over the real line^. Equiva- 
lently it may be written as 



F{x)=g[Ny{x)-N<{x) . 



(3) 



where Ny{x) {N^{x)) is the number of particles to the 
right (left) of x. The dynamics of this model, from vari- 
ous initial conditions and over different times scales, has 
been extensively explored in the literature (see references 
given above). 



C. Infinite system limit 

Let us consider now the infinite system limit, i.e., an 
infinite uniform distribution of points'^ on the real line 
with some mean density rig (^-g- a Poisson process). It 
is evident that the forces acting on particles are not well 
defined in this limit, as the difference between the num- 
ber of particles on the right and left of a given particle 
depends on how the limit is taken. Formally we can write 
the force field of Eq. (H) as 

F{x) = gno / dy sgn(y - x) + g / dy 5n{y) sgn(y - x) , 



(4) 

where 5n{y) = n{y) —no = J2i ^oiy — Xi) — uq represents 
the number density fluctuation. While the second term 
would, naively, be expected to converge if the fluctua- 
tions Sn{y) can decay sufficiently rapidly, the first term, 
due to the mean density, is explicitly badly defined (as 
the integral is only semi-convergent). Precisely the same 
problem arises for gravity in infinite 3 — c? distributions. 
The solution, known as the "Jeans swindle" , is the sub- 
traction of the contribution due to the mean density. As 
discussed by Kiessling in |Q , rather than a "swindle" , 
this is, in 3 — d, in fact a mathematically well-defined reg- 
ularisation of the physical problem, corresponding simply 
to the prescription that the force at a point be summed 
so that it vanishes in the limit of exact uniformity. The 
simplest form of such a prescription in 3 — d is that the 
force on a particle be calculated by summing symmet- 
rically about the particle (e.g. by summing about the 
considered point in spheres of radius R, and then send- 
ing R oo). This formulation needs no explicit use of 
a "background subtraction", since the term due to the 
mean density does not contribute when the sum is per- 
formed symmetrically. 



^ We use the standard convention that sgn(O) = 0, which impHes 
this same formula is vaUd for the force on a particle of the dis- 
tribution (rather than a test particle) at x. 

^ By "uniform" we mean that the point process has a well defined 
positive and spatially constant mean density, i.e., it becomes ho- 
mogeneous at sufficiently large scales. 



Applying the same reasoning to the 1 — d case would 
lead to the prescription 

F{x) = .9 J dy 5n{y) sgn{y - x) . (5) 

The question is whether this expression for the gravita- 
tional force is now well defined, and if it is, in what class 
of infinite point distributions. As we will detail in the 
next section of the paper, this question may be given a 
precise answer, as in 3 — d, by considering the proba- 
bility density function of the force in such distributions, 
described as stochastic point processes in infinite space. 
In the rest of this section we will simply explain the prob- 
lems which arise when the infinite system limit of expres- 
sion Eq. (^) is taken using a simple top-hat prescription. 
This discussion motivates the use of a smooth version of 
this prescription, which we then show rigorously in the 
subsequent section to give a well defined force for a broad 
class of infinite perturbed lattices. 

For Eq. (|^) to be well defined in an infinite point dis- 
tribution it must give the same answer no matter how it 
is calculated. Two evident top-hat prescriptions for its 
calculation are the following. On the one hand it may be 
written as 



F{x) 



g lim 

L — >oc 



x+L 



dy niy) sgn(y - x) 



(6) 



-L 



or, equivalently, 



F{x)=g lim N{x,x + L)- N{x- L,x) , (7) 

L — >oo L J 

where N[x,y) is the number of points between x and 
y, i.e., the force is proportional to the difference in the 
number of points on the right and left of x inside a sym- 
metric interval centred on x, when the size of the interval 
is taken to infinity. On the other hand, we can write 

r+L 

F{x) = g lim / dy Sn{y) sgn{y - x) , (8) 



L — ^oo 



or, equivalently, 

F{x) = g lim \n{x, L) - N{-L, x)] + 2g7iox, (9) 

i.e., we integrate the mass density fluctuations in a top- 
hat centred on some arbitrarily chosen origin. 

That these expressions are both badly defined in an in- 
finite Poisson distribution is easy to see: in this case the 
fluctuation in mass on the right of any point is uncorre- 
lated with that on the left, giving a typical force propor- 
tional to the square root of the mass in a randomly placed 
window of size L, which grows in proportion to ^/L (and 
thus diverges). Calculating the force with Eq. (^ one 
of us (AG) has shown in |g5[ that it is in fact not well 
deflned either in a class of more uniform distributions 
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FIG. 1: Calculation of the force using a top-hat regularisation 
centred on the point considered, i.e., as defined in Eq. (I). In 
an unperturbed lattice (case a) the force on points of the 
lattice vanishes. However, as shown in b) and c), when a 
single point is displaced off lattice, the force becomes badly 
defined, oscillating between g and zero as the size of top-hat 
goes to infinity. 



of points, randomly perturbed lattices''. Why this is so 
can be understood easily by considering, as illustrated in 
Fig. ^ the calculation of the force using Eq. (^ in such 
configurations. While on the unperturbed lattice (case a) 
the force on all points of the lattice is well-defined (and 
vanishing, as it should be), this is no longer true when a 
particle is displaced: the force on the displaced particle 
now oscillates deterministically (between g in case b, and 
zero in case c) and does not converge as L ^ oo. 

For the same case, of a single particle displaced off 
an infinite perfect lattice, the prescription Eq. (|^) for 
the force does, however, give a well-defined result if one 
chooses as origin a point of the unperturbed lattice: since 
the first ("particle") term is unchanged by the displace- 
ment of the particle, the only non- vanishing contribution 
comes from the second ("background") term, giving a 
finite force 



Fiu) = 2gnQU 



(10) 



where u is the displacement of the particle from its lat- 
tice site (i.e. the centre of symmetry) and we assume u is 
smaller than the lattice spacing. If we consider now, how- 
ever, applying random displacements of small amplitude 
(compared to the interparticle spacing) to the other par- 
ticles of the lattice, the problem of the first prescription 
Eq. (0) reappears: at any given L the first term in Eq. (||) 
picks up a stochastic fluctuation which varies discretely 
between ±g and zero, and does not converge as L ^ oo. 
This will evidently be the case for any such configura- 
tion generated by displacing particles off a lattice, and 



more generally for any stochastic particle distribution in 
1 — d. It is thus necessary to introduce some additional 
constraint to make this surface contribution to the force 
vanish. 

The previous literature on this model employ top-hat 
prescriptions equivalent to Eq. (||) to calculate the force, 
adding such a constraint. On the one hand, Aurell et 
al. in restrict themselves to the study of an infinite 
perfect lattice off which only a finite number are initially 
displaced. In this case the problematic surface fiuctua- 
tion vanishes for sufficiently large L. On the other hand 
]lO| , 16, |l^, |l^ impose exact symmetry in the displace- 
ments about some chosen point, which is then taken as 
the origin of the symmetric summation interval. A parti- 
cle entering (or leaving) at one extremity of the interval 
is then always compensated by one doing the same at the 
other extremity. 

We note that it is only in that the problem of the 
infinite system limit is actually considered. In the other 
works the authors do not discuss this limit explicitly: 
they consider and study in practice a finite system, with 
a prescription for the force equivalent to Eq. (Q) where 2L 
is the system size, i.e., without the explicit limit i — > oo. 
Symmetry about the origin is imposed because this al- 
lows one to use periodic boundary conditions. Such a 
finite periodic system of period 2L is equivalent to a fi- 
nite system of size L with reflecting boundary conditions. 
The dynamics of such a system is of course always well 
defined, for any (finite) initial distribution of the points 
in the box. This does not, however, mean that this dy- 
namics can be defined in the limit that the size of the 
system is taken to infinity. This is the question we focus 
on here, as the definition of such a limit is essential if a 
proper analogy is to be made with the cosmological prob- 
lem in 3 — d: in this case the gravitational force is well 
defined in the infinite system limit, for a class of statisti- 
cally translationally invariant distributions representing 
the initial conditions of cosmological models^ . 

The problems with the top-hat prescriptions arise, as 
we have seen, from non-convergent fluctuations at the 
surface of a top-hat window, which will be generic in sta- 
tistically translationally invariant point processes. It is 
thus natural to consider smoothing the summation win- 
dow, and specifically a prescription for Eq. (^ such as: 

F{x) = 5 lim / dy n{y) sgn{y - x) e"''!^"^! , (11) 



* The force is, however, shown to be well defined in this class of 
point distributions using the analogous definition for any power 
law interaction in which the pair force decays with separation. 
See H for details. 



Numerically one treats, of course, a periodic system, but it is an 
infinite periodic system, i.e., the force is calculated by summing 
over the particles in the finite box and all its (infinite) copies. 
This is the so-called "replica method" , used also widely in equi- 
librium systems such as the one component plasma teTJ. The 
infinite sum is usually calculated using the Ewald sum method. 
To obtain results independent of the chosen periodic box, the 
prescription for the force must converge in the appropriate class 
of infinite point distributions. 
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or, equivalently, 

F{x) = 5 lim y sgn(xi - a;)e-^l^'-^l , (12) 

i 

where the sum runs over aU particles in the (infinite) 
distribution. Rather than a smoothing of the summation 
window, this can be interpreted more physicaUy in terms 
of the screening of the gravitational interaction, i.e., the 
pair force law of Eq. (|l|) is replaced by 

/^(x) =-gsgn(x)e-''l-l, (13) 

and the gravitational force in the infinite system limit 
is defined as that obtained when the screening length is 
taken to infinity, after the infinite system is taken^ (see 
Fig. 2). This treatment is borrowed from the class of in- 
frared problems well known in quantum field theory. The 
standard procedure of handling infrared divergences is to 
apply an infrared regularization, to solve the regularized 
problem, and to remove the regularization at the end of 
the calculation, perhaps involving a renormalization. 





r 
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FIG. 2: Schematic representation of the smooth screening of 
the force (or, equivalently, summation window). 

For the case of a single particle displaced off a perfect 
lattice discussed above it is simple to calculate the force 
using Eq. (pi]). Denoting the lattice spacing by and 
the displacement by u, we have 

F{u) = 5 lim V sgn(n^ - M)e^^l"^-"l . (14) 

For |m| < ^ the sum gives 

2 sinh(/iu) ( ^ e-^"^) . (15) 

Expanding this in powers of /i we obtain 

^pH = ^ + 0(m). (16) 



Taking the limit ^ gives Eq. (|o|), i.e., the result ob- 
tained using the top-hat prescription Eq. (H) . The equiv- 
alence of the two prescriptions can likewise be shown to 
apply when displacements are applied to a finite num- 
ber of particles on the lattice (which leave the forces un- 
changed, and equal to Eq. (|l0|), if there are no crossings). 
Thus the only difference between the prescriptions is how 
they treat the contribution from particles at arbitrarily 
large distances when the infinite system limit is taken. 

We will show rigorously in the next section that, for a 
class of infinite perturbed lattices in which particles do 
not cross, the prescription Eq. (|l^) simply removes the 
problematic surface contribution present in the top-hat 
prescriptions (without applying any additional constraint 
of symmetry) . This gives a force on each particle equal to 
Eq. ( |lO| ) where u is the displacement of the particle, the 
only difference with respect to the case of a finite num- 
ber of displaced particles being that the origin of this 
displacement may be redefined by a net translation of 
the whole system induced by the infinite displacements. 
The force felt by each particle is thus equivalent to that 
exerted by an inverted harmonic oscillator about an (un- 
stable) equilibrium point. We note that this expression 
for the force is in fact what one would expect from a 
naive generalization of the analagous results in 3 — d. In 
the latter case it can be shown |2^ that the force on a 
single particle displaced off an infinite lattice by a vector 
u is, to linear order in |u|, simply 

F(u) = 47rG/9ou/3 . (17) 

This force is simply that which is inferred, by Gauss's 
law, as due to a uniform background of mass density -po 
(i.e. due to the mass of such a background contained in 
a sphere of radius |u|). The \ — d result is exactly anal- 
ogous, as 2no|u| is simply the mass inside the interval of 
"radius" While this result is valid, in 3 — d, only at 
linear order and for the case of a single displaced par- 
ticles, it is exactly valid in 1 — d in absence of particle 
crossings and for a broad class of displacement statistics. 
The reason is simply that in 1 — d the force on a particle 
is unaffected by displacements of other particles, unless 
the latter cross the considered particle. 

III. FORCES IN INFINITE PERTURBED 
LATTICES 

In this section we calculate, using the definition 
Eq. ( p^ , the gravitational force on particles in a class of 
infinite perturbed lattices. To do this we describe these 
point distributions as generated by a stochastic process 
in which the particles are displaced^. The force on a 



^ Although we will not use the interparticlo potential in our cal- 
culations, we note that /^(x) = —d4)p^/dx where (piii^) = 



-ge f^\^\/fM is the solution of ■ 



- n'^ij)^ = 2g5n{x). 



For an introduction to the formalism of stochastic point pro- 
cesses i.e. stochastic spatial distributions of point-particles with 
identical mass, see, e.g., p3]. 
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particle (or the force field at a point in space) is then it- 
self a stochastic variable, taking a different value in each 
realization of the point process, and the question of its 
defincdness can be cast in terms of the existence of the 
probability distribution function (PDF) of the force. We 
thus calculate here the PDF of the force on a particle 
with a given displacement u, in the ensemble of realiza- 
tions of the displacements of the other particles. The 
result is that, for the class of stochastic displacement 
fields in which displacements are such that particles do 
not cross, this force PDF becomes simply a Dirac delta 
function. This gives the anticipated result, that the only 
force which results is that due to the particle's own dis- 
placement given by Eq. (^0|), modulo an additional term 
describing a contribution from the coherent displacement 
of the whole infinite lattice if the average displacement is 



non-zero. 



A. Stochastic perturbed lattices 

Let us consider first an infinite 1 — d regular chain of 
unitary mass particles with lattice spacing £ > 0, i.e., the 
position of the n*'' particle is X„ = n£, and the micro- 
scopic number density can be written as 



l^in (^) 



+ C30 



(5_D (x — n£) . 



(18) 



We now apply a stochastic displacement field {[/„} to 
this system, in which the displacement J7„ is applied to 
the generic n*'' particle with n G Z. Let us call 
the single realization of the stochastic field {[/„}. The 
corresponding realization of the point process thus has 
microscopic number density 



n{x) 



+00 

E 

71 — — OO 



(19) 



This displacement field is completely characterized by the 
joint displacement PDF 'P{{un}) where is the set 

of all particle displacements with n e Z. We will fur- 
ther assume that this stochastic process is statistically 
translationally invariant, i.e. V{{un}) — V{{un+i}) for 
any integer I. This implies in particular that the one dis- 
placement PDF (for the displacement applied to a single 
particle) is independent of the position of that particle, 
i.e., the function 



B. Mean value and variance of the total force 

Let us denote in general by Ffj^{xo) the total gravita- 
tional force, with finite screening fi, acting on the particle 
at xq and due to all the other particles placed at x„: 



F^L{xo) = .9^sgn(a 



xo)e 



'fj.\x„-xo\ 



(21) 



Writing now Xn = ni + u„ in Eq. ( pT] ) , we can write the 
total screened force on the particle at xq = uq in a per- 
turbed lattice for a given realization of the displacement 
field: 

F^{uo)=gY,sS^i^^ + ^n~uo)e-^\"'+""-^"\. (22) 

Note that, given the assumed statistical translational 
invariance of the field {Un} the statistical properties 
of the force are the same for all particles in the sys- 
tem. If, further, we assume now that the displacements 
from the lattice are such that particles do not cross, i.e. 
sgn(ri^-|- w„ — uq) — sgn(n) for n ^ 0, this can be written 
as 



(23) 



where we define for, n > 1, 
fn = /«(m) 



g-/i(n„-Mo) _ g-/i(-UO-"-n) 



We now take the average of Eq. (^3|) over all realiza- 
tions of the displacements of all particles, except the cho- 
sen one uo, which we consider as fixed. We denote this 
conditional average as (•)q, while we use (•) for the uncon- 
ditional average. In order to do this we need the condi- 
tional PDF of Un to Uq, which by definition of conditional 
probability is 



Pn(u;uo) 



qnQ[U,Uo) 



(24) 



Piuo) 

By using this function we can write 

(/n(M))o = e''^°Pni^l;uo) - e-^"«P_„(-/x; wo) (25) 
and therefore 



[[dunV{{un})SDiu - w„) (20) (F^iuo))^ = 5^ U"'"' Pn{f^;uo) - e^^"«P_„(-^; mq) 



is independent of m, i.e. Pm{u) — p{u). Moreover the 
joint two-displacement PDF 

qnm{u,v) = / '^dUn'P{{Un})5D{u - Um)6D{v - U„) 
•' n 

depends parametrically on the lattice positions n, m only 
through their relative distance (m — n) . 



n=l 

where we have defined 



-Pn(M;wo) 



E 

k=0 



duPn{u;uQ)e 

ki 



(26) 
(27) 
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The latter equality is valid when all the moments (C/,i)q 
of Pn{u;uo) are finite. Note that, given the assumption 
that particles do not cross, it follows from the definition 
( p^ ) that 5„o(w,wo) — for u + n ^ uq respectively for 
n ^ 0. Therefore Pn{u;uQ) is always zero for some suf- 
ficiently negative uq dependent value of u if n > 0, and 
likewise for sufficiently positive values if ti < 0. This 
ensures that the integral in Eq. ( p7|) is indeed finite. 

In order to study the behavior of Eq. ( p6| ) for /i — ^ 0, 
we will assume that 



\n—m\ — ^oo 



p{u)p{v) . 



(28) 



This corresponds to the assumption that the displace- 
ment field is a well defined stochastic field, which requires 
(see e.g. [^) that the two-displacement correlations van- 
ish as the spatial separation diverges. We will discuss in 
the next section the restriction this corresponds to on the 
large scale behaviour of the density perturbations, which 
is of particular relevance when one considers the analogy 
to 3 — d cosmological simulations. 
Assuming Eq. (ESh we can write 



Pn{u; uo) = p{u) + r„(u; uq) , 

where r„(M; uq) is a function vanishing for |n| oo and 
with zero integral over u for any n. As a consequence 



P„(/i; uo) = p{^^) + fnifJ.; uq) 



(29) 



where we used the definition analogous to Eq. ( p7| ) for 
p(/i) and f„(/i;Mo), and the latter vanishes for jj. 
and/or n — *■ oo. If we now suppose that both (U) and 
(C/„)g are finite, with evidently (J7„)q (U) for n oo, 
we can write at lower order: 

pi^l) = ^-^i{u) + o{^l), (30) 

f„(/X; Uq) = /i((C/) - (C/„)o) + o(/i) • 

It is now simple, by substituting Eqs. (|9|) and ( |30| ) into 
Eq. (|2^), to show that, if ( (U) — (C/„)g) decays in n as 
a negative power law or faster, we have 

{F{uo)), = lim (F^iuo)), = 2gno{uo - (U)) . (31) 

We will now show that both for uncorrelated displace- 
ments, and then more generally for correlated displace- 
ments with decaying correlations, this average force is in 
fact the exact force in every realization. We do so by 
simply showing that 



lim 



0. 



(32) 



This implies that the variance of the conditional PDF of 
the total force F acting on the particle at uq vanishes, 
i.e., this PDF is a Dirac delta function at the average 
value given by Eq. (^). Compared to the simple case of 
a single displaced particle we analysed above, the only 
effect of the (infinite number of) other displacements is 



to possibly shift the centre of mass of the whole (infinite) 
distribution with respect to which the displacement of 
the single particle is defined. 

In order to show Eq. (|3^ ) wc note first that the second 
conditional moment of F may be written 

1,C30 

n,m 

oo 

= {F,{uo))l+9'Y.e-'^"'A,,{f,) 

l,oo 

(33) 



with 



M^^) = {fn)o-{fn)l 
BnmifJ') = (/n/m)o - (/«)o (/m)o 7^ n) , 



(34) 



and where ^ as usual indicates the sum over m and ; 
with the exception of the n ~ m terms. To prove Eq. (32 
it is sufficient to show that the last two terms in Eq. (33 
go continuously to zero as /i does so. 



C. Lattice with uncorrelated displacements 

We consider first the case that the displacements are 
uncorrelated and identically distributed, i.e.. 



-t-oo 



(35) 



n— — oo 



We refer to this as a "shuffled lattice" configuration (fol- 
lowing (22I). In this case conditional and unconditional 
averages coincide. Given the assumption that the dis- 
placements do not make particles cross, we must have 
that p{u) — for \u\ > £/2, implying that all the mo- 
ments of p{u) are necessarily finite. 

In this case the w„ are statistically independent and 
identically distributed random variables. Given the defi- 
nition Eq. (^4|), it follows that the /„ also have this prop- 
erty, i.e.. 



(/n/m) — ifn) (/m) 



(36) 



and thus that i?,im(/i) = 0. Further An{fi) is independent 
of n and can be expressed explicitly as 

An(p) = e'^-° [p{2^^yp'{^^)] -e^^A-^o [p(_2^)_p2(_^)] 

(37) 

Expanding this expression in fi about /i = 0, we find 
that the leading non- vanishing term is at order /i^. The 
desired result, Eq. (|3^), follows as 



n=l 



where 0{^^) means as usual a term of order I in /i. 
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D. Lattice with correlated displacements 

We now consider the case where the displacements are 
non-trivially correlated. In order to calculate A„(/z) and 
BnmifJ') we need both the conditional single displacement 
PDF Pn{u;uo) and the conditional two-displacement 
PDF (5„m('u, w; Uq), both conditioned to the fixed value 
Mo of the stochastic displacement Uq. The function 
Qnm{u, v; Uq) is defined by the rules of conditional prob- 
ability as 



Snmo{u,V,UQ) 

p{uo) 



where Snmi(u,v,w) is the joint three displacement PDF 
of having the three displacements u, v, w respectively at 
the lattice sites n,m,l. 

Let us start from the evaluation of A„(/i). From its 
definition it is simple to show that 

-2Q„_„(^,-^;uo), (38) 



where 



+00 



In order to study the limit n —^ we have to expand 
Uq) and Qmyiip-, ±/^; uo) in powers of ^. Assuming 
that at least the first two moments of the displacement 
statistics are finite, we can write 

+ (C/2„)^±(C/„C/,„)o) 

(39) 

Using this result and Eqs. ( |25| ) and ( ^ ) in the definition 
(^) of A„(/i), it is simple to show that 



(40) 



+2 {{UnU^n)^ - (CZ„)o (C^-n>o)] + ■ 

Note that for |n| — > cx) we have (t/„)o — > (J^), (^^,1)0 ~* 
(C/2) and (C/„C/_„)o (C/)^ Therefore we can write 

An{^x) m'(((:^') - ((7)')(e2^"« -I- e-2M«o) ^ 

where we have used the fact that, as the coefficients of 
the higher order contributions in /i to A„(/i) are non- 
diverging, they can be neglected. This is sufficient to 
conclude that 



Y,e->^■^A^{^i)^o{^i), 



(41) 



where 0{fJ) as usual means a term of order /x', and there- 
fore the sum vanishes as fi for /i — > 0. 



Let us now move to analyze the last sum in Eq. 
We study the behavior of -B„m(/i) as defined by Eq. 
It is simple to show that 



H). 



(/n/m)o = e 2^''«Q„„(/x,/i;uo) + e2^"''Q-«-m(-M,-M;wo) 

-Qn--m(M7-M;'"o) - Q-„m(-M7/^;Uo)- (42) 

Using this equation together with Eqs. (^,(|25|) and (|3S|), 
we can write 

BnM = M'[e"'''""g(","^;^io)+e'^"°.g(-n,-m;iio) 
-5(n, -m; mq) - 5(-"-: "o)] + o(a<^), 

(43) 

where we have called 

g{n,m;uo) = (C/„C/„i)o - (f^n>o (f^™>o ' 

i.e., the conditional displacement covariance matrix. 
Since this is a "conditional" correlation it does not de- 
pend simply on n — TO, but on both n and to in a non- 
trivial way. However for both |n|,|TO| — )■ oo the condi- 
tional averages coincide with the unconditional ones and 
therefore we can write 



g{n, to; uq) = c{\n — m\)[l + h{n, to; uq)] 



(44) 



where c(|n — to|) = (UnUm) — (U)'^ is the unconditional 
displacement covariance matrix, and h{n,m;uo) — > for 
|n|, |to| oo. In order to analyze the asymptotic behav- 
ior for small fi of 



1,00 



/(/i)EE^'e-''("+")B„,„(Ai), 



(45) 



it is sufficient to study the behavior of the sum coming 
from the first term (or equivalently the second) of Bnm (m) 
in Eq. (^) as it is the most slowly convergent one, i.e., 
basically to study the following sum: 

1,C30 

J(/.) = ^'e-^("+'")5(n,m;zio). 



Since h{n,m;uo) for |n|, \m\ 00, the small /j, scal- 
ing behavior of J(/i) is the same if we replace g{n, to; uq) 
by c{\n — to|): 



j(m) 



1,00 

n.m 



(|n — m|) . 



(46) 



This can be also shown by the following argument: as- 
suming that /i(n, to;uo) is bounded, say |/i(rt, to; uo)| l£ 
A, we can write 

l-/(M)l<El:™'e-^("+'")|5(ri,m;«o)| 
<(l + ^)Ei™'e-^^"+"^|c(l?^-m|)|. 
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Therefore the convergence to zero of fi^ times the right- 
hand side of Eq. ( ^ ) is a sufficient condition to have the 
variance of F to vanish for /i 0. 

Let us now analyze the right-hand side of Eq. (|4^) . We 
can write 

1,00 

^'e-''("+'")c(|n-m|) 

1,CXD 

^^,-M"+™)e(|n-m|)-c(0)^5^, (47) 



i.e. it vanishes in the ^ hmit and the PDF of the 
total force acting on a particle displaced by uq from its 
lattice position is W{F;uo) = S[F ~ 2g{uo — (C/))]. In 
other words, even in the case of spatially correlated dis- 
placements, the total force acting on a particle is a de- 
terministic quantity equal to 2g{uQ — (U)) with no fluc- 
tuations. This value depends only on the displacement 
of the particle on which we are calculating the force and 
not on the displacements of other particles as it does in 
3-d |26l. 



where c(0) is the single displacement variance. Note that 
the second term is of order at small ^ and therefore 
gives rise to a term at linear order in /i in Eq. (^). Let 
us introduce the Fourier transform c{k) of c(n), defined 

by 



dk 



iky 



Using this in the right-hand side of Eq. (E^ we get 



^e-^("+™)c(|n-m|) 



(48) 



_^ 27r ^ e^A" + 1 - 2et^ cos k 



The small /i limit of this integral is dominated by the 
behavior at small k of the integrand. In this limit the 
following approximation holds (e^^ + 1 — 2e^ cos k) ~ 



(^ -I- A: ). Let us also assume that c{n) 



at large n 



(with in general a > 0) which implies at small c{k) ~ 
|fc|"~^ for < a < 1 (with logarithmic corrections for 
a = 1) and c(k) ^ \k\^ with (3 > for a > 1. Therefore 
the small fi behavior of Eq. (^) is the same as that of 
the simple integral 



' dk c{k) 

, 2^/i2-Hfc2 



/i'^^^ for a > 1 



for < a < 1, 



(49) 



Taking also into account the second term in Eq. 
can therefore conclude that 



^'i?„™(A*)e'(" 



for < a < 1, 
/i for a > 1 . 



(50) 



This, together with the results for the first sum in 
Eq. (^), it follows that at small /i 



fi" for < a < 1, , ^ 
fj, for a > 1, ^ ' 



The case of a decay faster than any power, e.g. exponential 
decay, can be included for a — > 00. 



IV. DYNAMICS OF ID GRAVITATIONAL 
SYSTEMS 

In the previous section we have shown the prescription 
Eq. (^ij) for the 1 — d gravitational force to give a well 
defined result in a class of infinite displaced lattice dis- 
tributions. This result can be used in the construction 
of different toy models, through different prescriptions 
for the dynamics associated to these forces. In this sec- 
tion we discuss two such models, analogous to the 3 — d 
cases of gravitational clustering in an infinite static or ex- 
panding universe, respectively. In the last subsection we 
discuss in detail the relation of these models to previous 
treatments of such models in the literature. 

As motivation let us first comment on the reason for 
our interest in the case of perturbed lattices: in 3 — d 
cosmological A'^-body simulations precisely such configu- 
rations are used as initial conditions. The reason is that 
by displacing particles from a lattice in this way, one 
can represent accurately, at sufficiently large scales, low- 
amplitude density perturbations about uniformity with 
a desired power spectrum P{k) (for a detailed discussion 
see e.g. |2^ or ^^l)- This algorithm is strictly valid in 
the limit of very small relative displacements of particles, 
so that the assumption that particles do not cross in our 
derivation is a reasonable one (although not, as we will 
discuss in our conclusions, rigorously valid). The fur- 
ther assumption Eq. (|2|) we have made, on the decay of 
correlations, corresponds, also to a reasonable restriction 
on the class of initial power spectra. Indeed it can be 
shown easily that it corresponds, in d dimensions, to the 
assumption that P{k)/k'^ be integrable at fc = 0. In 3 — c? 
this corresponds to P{k ^ 0) ~ A:" with n > — 1, which 
is strictly satisfied in typical cosmological models which 
are characterised by an exponent n = 1 at asymptotically 
small fc. 



A. Toy models: static 

The simplest such model is the conservative Newtonian 
dynamics associated to the derived force law, i.e., with 
equation of motion 



Xi = F^{{xj,j = O..oo},i), 



(52) 
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where Fi is the gravitational force on the i-th particle 
of the distribution, with position Xi at time t (and dots 
denote derivatives with respect to t), calculated using the 
prescription Eq. (|l^, i.e., 



static universe, with equations of motion 



-g lini ^ sgn{xi - Xj)e ^1^' ' . 



(53) 



3^1 



We have shown that, for the case of an infinite lattice 
subjected to displacements which (i) do not make the 
particles cross, and (ii) satisfy Eq. (|2^), the force on the 
right-hand side is simply given deterministically as pro- 
portional to the particle's displacement (when (C/), the 
average displacement, is zero). Denoting then the dis- 
placements of the z-th particle by Ui, i.e. Xi — ii + Ui, 
the equation of motion is therefore 



(54) 



i.e., simply that of an inverted harmonic oscillator. The 
same equation is valid in the case that (U) ^ if we de- 
fine Xi = i£+ (U) +Ui. This equation of motion is valid, of 
course, only as long as the non-crossing condition is satis- 
fied. While it is in principle straightforward to generalize 
our calculation of the force to incorporate the effects of a 
finite number of crossings, it is much more convenient to 
make use of the following fact, which we recalled above: 
particles crossings in 1 — d are equivalent, up to exchange 
of particle labels, to elastic collisions between particles, in 
which velocities are exchanged. This means that if we are 
interested in properties of the model which do not depend 
on particle labels, the model of 1 — c? self-gravitating par- 
ticles is equivalent to a model in which particles bounce 
elastically. In this case the particles displacements from 
their original lattice sites are at all times such that there 
is no crossing of particles, and Eq. ( p4[ ) remains valid, ex- 
cept exactly at "collisions" . The dynamics of this model 
is therefore equivalent to that of an infinite set of inverted 
harmonic oscillators centred on the sites of a perfect lat- 
tice which bounce elastically, exchanging velocities, when 
they collide. As in the finite "sheet model" the equation 
of motion may be integrated exactly. Defining, for con- 
venience, time in units of the characteristic "dynamical" 
time Tdyn = 1/ V2gno J the evolution between collisions is 
given exactly by 

Ui{tf) +t)~ Ui{to) coshi -I- Vi{to) sinhi, (55) 
Vi{tQ + t) — Uiito) smht + Vi{tQ) cosht, (56) 

where Ui(to) [vi{to)] is the position (velocity) after the 
preceeding collision. The solution of the dynamics re- 
quires simply the determination of the next crossing time, 
which involves the solution of a quadratic equation (in 
e*), followed by an appropriate updating of the velocities 
of the colliding particles. 



r,; - r,, 



(57) 



for identical particles of mass m. We use the superscript 
J on the sum to indicate that the sum is calculated using 
the Jeans swindle. As we have discussed this "swindle" 
in 3 — d can be implemented by summing symmetrically 
about the point i either in a top-hat (i.e. sphere) or using 
the limiting procedure with a screening. 

The equations of motion for particles in an infinite ex- 
panding 3 — (i universe are usually written in the form 



2Hx, = -- 



.J X,- 



.3 ^ 



(58) 



where are the so-called comoving coordinates of the 
particles, H{t) = a/a is the Hubble "constant", and a{t) 
is the scale factor which is a solution of the equations 



(59) 
(60) 







SnG 
3a3 ^ 


C 


d 
a 


AttG 
3a3 







where po is the mean mass density when a = 1, and C is 
a constant of integration^. 

Note that these equations can be derived entirely in a 
Newtonian framework, and correspond simply to a differ- 
ent regularisation of the infinite system limit than that 
employed in the Jeans swindle: instead of discarding the 
effect of the mean mass density, the force is regularised so 
that the mean density sources a homologous expansion 
(or contraction) of the whole system. This corresponds 
to taking equations of motion 



-Gm lim - — ^, 

B^oo ^ r, — r,- ' 



(61) 



i.e. with the sum for the force calculated by summing 
symmetrically about a chosen origin. Dividing the sum 
into a term due to the mean mass density and a term due 
to fluctuations about this density, this may be written as 



AnGp 



GmY^ 



J r, 



(62) 



Neglecting the second term (i.e. taking only the force 
due to the mean density) gives an equation of motion 
admitting solutions of the form ri{t) = a(i)rj(io), with 



B. Toy models: expanding 

The model we have just discussed is the 1 — d analogy 
for the problem of gravitational clustering in an infinite 



^ C = corresponds to the fiat Einstein de Sitter universe, C > 
to a closed universe, and C < to an open universe. In the 
Newtonian derivation of these equations, given below, C can be 
expressed in terms of the physical particle velocities at some 
initial time. 
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a{t) satisfying Eqs. ( |59| ) and (pOj). Changing to comov- 
ing coordinates defined by r,j = a{t)xi in Eq. (|6l| ) [or in 
Eq. (|6^)], and using Eq. (|60|), then gives Eq. (|58|). 

Note that setting a(t) = 1 in Eq. ( jSq ) gives exactly the 
static case Eq. ([S^), i.e., the Jeans swindle in static space 
corresponds formally to the non-expanding limit of an ex- 
panding FRW universe. This static solution a{t) — 1 is, 
however, a solution to Eqs. (59) and ( |6^ ) only ii po — 
(and C = 0), i.e., it is not a physical limit of the ex- 
panding case but corresponds to the different prescrip- 
tion, Eq. (57), for calculating the force in the infinite 
volume limit. While almost all numerical studies are of 
the expanding case (for a review, see e.g., [Q), a recent 
study of the static case for such initial conditions has 
shown that the evolution of clustering is, in essential re- 
spects, qualitatively similar in both cases. This suggests 
that it may be possible to understand essential qualita- 
tive features of the dynamics of structure formation in 
the universe in the conceptually simpler framework in 
which there is no expansion. 

With the 3 — d equation of motion in the form of 
Eq. (|5^) it is evident how the static 1 — d model dis- 
cussed above is naturally modified to mimic the 3 — d 
expanding case: one can simply replace the force term 
due to the infinite 3 — d distribution [i.e. the sum on the 
right-hand side of Eqs. (^8|)] by that due to the 3 — d 
distribution consisting of infinite sheets. The summation 
prescription implementing the Jeans swindle for the gen- 
eral 3 — d case, i.e. spherical top-hat summation, is then, 
as we have discussed at length above, most appropriately 
replaced by the smooth prescription we have given. Thus 
we take the following 1 — d equation for the positions Xi 
of the particles (sheets): 



2Hx, 



27rGI] 



lim > sgn(a 



(63) 

where the sum extends over the infinite distribution of 
sheets, and we have explicitly made the identification g = 
27rGI] (where E is the mass per unit surface). 

With initial conditions in the class of 1 — d infinite 
perturbed lattices for which we have shown the sum for 
the force to be well defined and given by Eq. (|o|), we 
then have 



2Hui 



AttGpq 



(64) 



where we have used that the mean comoving mass density 
Po = Eno (i.e. physical mass density when a = 1). As in 
the static case, this equation of motion remains valid at 
all times if we exchange the labels of particles when they 
cross, so that they bounce instead of passing through one 
another. 

For the case of an Einstein de Sitter (EdS) uni- 
verse, which corresponds to C = in Eq. ([59[), a{t) — 
(67rGpo)^/^i^/^ and Eqs. iM) simplify to 



3t 



3*2 



(65) 



of which the independent solutions are Ui{t) oc t^/s and 
Ui{t) (X [which are simply the well known growing 
and decaying solutions for small perturbations to a self- 



gravitating fluid in an EdS universe (see, e.g., [^ )] 
evolution in between "collisions" is thus given by 



The 



Ui{t) = Uiito) 



+V^{to)ta 



2/3 



2/3 



(66) 



Note that, from Eq. (|66|) the determination of the cross- 
ings in these models, instead of a quadratic equation in 
the static model, thus involves the solution of a fifth order 
equation (for t^/"^). 



C. Discussion of previous literature 

1. Static models 

A few previous studies [|l^, |l^, |l^ have considered 
static 1 — d toy models, defining the force on the right 
hand side of Eq. (^2|) as the derivative of a potential, 
which is the sum of the contribution from the sheets in 
a finite system of size L, and an additional one due to a 
uniform negative background. This is exactly the "naive" 
version of the Jeans swindle discussed above, and corre- 
sponds exactly to the prescription Eq. (||) for the calula- 
tion of the force (with L finite) . The authors of [^ dis- 
cuss explicitly the problems associated with taking the 
infinite system limit. As a result they limit their analysis 
only to a case for which their prescription gives a unique 
and finite result: a finite number of particles displaced 
off an infinite perfect lattice, modelling a finite localized 
perturbation embedded in an otherwise uniform universe. 
It is simple to verify that equation of motion for these 
displacements is then exactly Eq. (54), which we have 
now shown to be valid for the infinite lattice with per- 
turbations which do not break the lattice translational 
invariance. 

In [|l^, , on the other hand, the dynamics is formu- 
lated for a system of finite L, and the problem of the 
definedness of the force in the infinite system limit is not 
explicitly addressed. Instead it is dealt with implicitly 
by assuming that the finite system is symmetric about 
some point. Taking this latter point as origin of coor- 
dinates, the top-hat prescription Eq. (^ for the force at 
coordinate position x may then be rewritten as 



F{x) = ~2gN{0, x) + 2gnox , 



(67) 



in which the size of the system does not explicitly appear. 
Labelling the particles by their position with respect to 
the origin {i — 1...N), the force on the i-th particle may 
then be written 



Fi = 2gno 



- ( ^ )(* 



1 



(68) 
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FIG. 3: The variance of the displacement Ui (see text) as 
a function of a particle's ordered position i, calculated for 
one thousand realizations of one thousand particles randomly 
placed in an interval. 



a function of i, as measured in one thousand realizations 
of one thousand randomly thrown particles. In a typical 
realization the force on a particle in the centre of the box 
is thus much larger than that on a particle at the bound- 
aries. This means that the typical force on a particle 
therefore not only diverges as L does, but that in a finite 
system its typical value depends on the position of the 
particle with respect to the boundaries. In practice this 
means that all the particle trajectories in a symmetric fi- 
nite system of initially Poisson distributed particles are, 
right from the initial time, modified when L is increased, 
and do not converge well as L — *■ oo. Simulations of such 
initial conditions reported in |l^ show the associated co- 
herent global evolution of the mass distribution, which 
contrasts qualitatively with the local clustering charac- 
teristic of the I — d (and cosmological) simulations which 
we will describe in the next section. 



2. Expanding models 



where Xi is the position of the particle. For any finite 
system the quantity in brackets can be considered as the 
displacement Ui of the particle i from its "original" lattice 
site [at {i — 1)L/N]. Thus the equation of motion for the 
particles is again identical to that we have derived. 

We note again that we have derived this force law in 
this article without the assumption of symmetry (and 
without the explicit introduction of a background). Fur- 
ther, and most crucially, we have shown it to remain valid 
for a certain class of distributions when the infinite vol- 
ume limit is taken — perturbed lattices without crossing 
and displacements of finite variance. In this respect we 
underline, as we have done in Sect. ||, that while in the 
formulation of Q the same equations of motion Eq. ( |5^ ) 
are valid for the particles in any finite symmetric system, 
this does not mean that the infinite system limit is well 
defined, even with the assumed symmetry. It is illustra- 
tive to see what "goes wrong" when the infinite system 
limit is taken, specifically, for the case of a Poisson dis- 
tribution, i.e., when we consider a system of size L in 
which we distribute N particles randomly, and then take 
L — > oo at fixed uq = N/L. The problem is that forces, 
although defined at any finite L by Eq. (|68|), diverge as 
L does. This can most easily be seen by considering the 
force written as Eq. (|67|): the force on a particle at x, 
as it is proportional to the fluctuation in the number 
of particles in the interval [0, x] about its average value, 
grows in proportion to ^/x. Working with Eq. (68) this 
corresponds to the fact, which can easily be shown, that 
the variance of the displacements Ui (as defined above) 
diverges, violating an essential assumption for the per- 
turbed lattices in the preceeding section. Further this 
variance (in the finite system) depends also on «, so that 
discrete translational invariance, which we also assumed, 
is likewise broken. These properties are illustrated in 
Fig. 1^, which shows the variance of the displacement Ui as 



We note first that Eq. ( |64|) coincides exactly with that 
obtained in the so-called Zeldovich approximation in 3— d 
(see, e.g., |2^), when Ui is replaced by a vector func- 
tion u(x). This approximation describes the evolution 
of displacement fields u(x) engendering small amplitude 
fiuctuations to a self-gravitating fluid in an expanding 
universe, and can be obtained rigorously by a pcrtur- 
bative treatment of the full fluid equations ||2^ in the 
lagrangian formalism^°. For the case of one-dimensional 
perturbations it is well known (see | pT| and references 
therein) that this approximation becomes exact, up to 
the time when caustics form, corresponding to the cross- 
ing of "sheets" of fluid (i.e. particles in our case). It 
is thus, perhaps, not surprising, a posteriori, that we 
recover exactly the Zeldovich approximation for the mo- 
tion of discrete sheets up to the time they cross: as the 
pair force between sheets is independent of separation, 
the only way a sheet can "see" that the force sourcing 
its motion is discrete, rather than continuous (as in the 
fluid limit), is when it crosses other sheets. 

Eq. (64) can equally be derived [E^ Is^] using a pertur- 



bative treatment of the dynamics of an inflnite perturbed 
lattice (in 3 — d) of particles. For plane wave displace- 
ments of the particles with a wave- vector orthogonal to 
one of the lattice planes, the amplitude of the displace- 
ment wave obeys exactly this equation in the limit that 
the discreteness of the mass distribution in these orthogo- 
nal planes is neglected. This latter assumption is weaker 
than that used in this framework to derive the Zeldovich 
approximation for a general perturbation, which would 
require also that the displacement be of long wavelength 
compared to the discreteness scale in the direction par- 



X is a lagrangian coordinate and the fluid is exactly uniform when 
u(x) = 0. 
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the authors study exactly 
( |64|) for the displacements 



allel to it. 

In the studies of 
the equations of motion Eq. 
of sheets perturbed off a perfect lattice (as in cosmolog- 
ical simulations). They adopt these equations arguing 
that they represent the fluid limit for 1 — d perturbations 
in a 3 — d expanding universe. While before sheet crossing 
(i.e. the formation of caustics), as discussed above, this 
is indeed known to be true — these equations are just 
the Zeldovich approximation which is, in this regime, ex- 
act — the extension to longer times is argued to be valid 
because the "collisionless" sheets of fluid will simply pass 
through one another. Our derivation of these equations 
shows that this in fact corresponds to the discrete par- 
ticle/sheet model. Indeed we have not taken the fluid 
limit in our derivation, and the equations do not repre- 
sent the fluid limit of this model. It simply happens to 
be the case that in this model, before crossing, the equa- 
tions correspond with those in the fluid limit, for the 
physical reasons we have mentioned above. After cross- 
ing this equivalence breaks down, and the prescription 
used by |1^ to "analytically continue" the fluid model 
beyond its regime of validity actually maps onto this dis- 
crete particle/shell model. We will discuss further in our 
conclusions the definition of the fluid limit for the I — d 
models we are considering here. 

The other two groups who have considered 1 — d toy 
models incorporating 3 — c? expansion have, as in this 
article, worked in a particle/sheet framework. Both the 
original model, proposed in ||l0|| and studied further in 
1 14 1, and the subsequent one proposed and studied in 
1 15, |l^ and |l^, derive their (different) equations of 
motion by following, formally, the steps described above 
leading from Eq. (|6l|) to Eq. (pq). The force on the right 



hand side of Eq. (57) is simply that due to the sheets, 
calculated in the analogous manner^^, i.e.. 



fi ~ 2tiGT. lim 

L- 



^ sgn(r 



(69) 



The change to comoving coordinates, when assumed also 
to rescale the mass in the sheets in the orthogonal direc- 
tion (so that S — + S/a^), gives 



lim y sgn{xj — Xi) + 2nQXi 

Xje[-L.L] 



(70) 



In [0 the force term is simply denoted Ei, without an expUcit 
prescription for calculating it. It can be inferred from the de- 
scription given subsequently of the numerical simulations that 
the implicit summatiori is the symmetric top-hat centred at the 
spatial origin. In [ |L5| , |16[ , on the other hand, the top-hat regu- 
larisation is explicited. 



provided that a{t) obeys the equation 

a _ 47rG 
a 



-Po ■ 



(71) 
d) when 



As above po — Sno is the mass density (in 3 
a = 1. 

The Eqs. (|7^) are those adopted by 0, |l|, |l|, |l|, 0. 
The term which we have written on the right hand side 
of the equation corresponds exactly to the prescription 
Eq. (^ for the calculation of the force. It incorporates 
the required subtraction of the effect of the background, 
so that motion in comoving coordinates is sourced only by 
perturbations to uniformity. Just as in the static models 
of 1 13, |l8| discussed above, which are obtained formally 
by setting a = 1 in Eq. ([TOj), this force is well defined 
only if symmetry is assumed about the chosen origin in 
the point distribution. This is indeed the assumption 
made in the numerical studies of jlO, ^ p7[ . 

The difference between the mode s of ]llO|, |14|| and of 
p^ , (studied also in |l^) arises only in what is as- 
sumed in each case about the scale factor a{t). The for- 
mer authors impose an EdS cosmology behaviour for the 
scale factor, a{t) cx i^/'^, and require that it is a solution 
of Eq. (|7l|). While this is mathematically consistent, it is 
not physically coherent: comparing Eq. ( |7l] ) and Eq. ( [59| ) 
we see that it corresponds to imposing a Hubble expan- 
sion sourced by a mean density three times the physi- 
cal mass density of the sheet (or, equivalently, assuming 
that the gravitational constant is not the same for the 
background as for the perturbations). Refs. ||l5|, |l^, on 
the other hand, simply impose that a(t) be the EdS ex- 
pansion, with the right normalization. This amounts to 
adding "by hand" a term to the derived equation [ p7[ . 
While this is not mathematically rigorous within the con- 
text of the derivation just described, it is more appropri- 
ate physically. Indeed it corresponds effectively to simply 
replacing the Jeans swindle 3 — c? force term in Eq. ( |5§| ) 
by the prescription Eq. (^ . This differs from the "deriva- 
tion" we have given above for Eqs. ( |63| ) only in the form of 
the Jeans swindle adopted. For the case that symmetry 
about the origin is assumed, we have the same equations 
of motion. In a finite system Eq. ( |68| ) is valid and so 
the equations of motion in their numerical simulations 
reduce exactly to Eqs. (p^). 

In conclusion the equations of motion Eqs. ( |6^ ) are 
exactly the same as those used by 0, 0, and by 
]l5| , [iq , [iTt . The only difference in practice between all 
these studies are the initial conditions adopted and also 
the analysis of the resultant clustering given. Rather 
than working in the cosmological time variables, the lat- 
ter authors define, a new time coordinate r = •y/2/31nt. 
Eqs. (|65|), for the case of an EdS universe, then take the 
very simple form 



dr2 



1 dui 
7^ 



(72) 



In these variables the model is thus equivalent to an in- 
finite set of inverted oscillators which bounce elastically. 
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with an additional constant damping. Because of the 
fifth order equation which must be solved to determine 
the crossings (now for the parameter t^^^ = e"^/^), the 
model has been dubbed the "quintic" model by the au- 
thors of @. 

The model of ||l^, on the other hand, is mathe- 
matically consistent but of less apparent relevance to the 
"real" 3 — d cosmological model, as the expansion it im- 
poses on the sheets is not the physical 3 — d expansion. 
Indeed in this respect we note that, in the derivation of 
[0, any function a{t) satisfying Eq. { fn\ ) can be adopted 
with the same consistency. The only way in fact in which 
the derivation of the 3 — d equations can be rigorously 
adapted to 1 — d is by using the 1 — d expansion law de- 
rived from Eq. (^9|) in the limit of uniformly distributed 
sheets. This is 



a{t) = 1 + Hat- 27rGnot^ 



(73) 



where Hq = H{t — 0), i.e., free fall in a constant gravi- 
tational field of strength AiiGnQ. As this is very different 
to the 3 — d expansion law it is probably not a variant of 
the toy model which is of practical interest. 



V. CASE STUDY: INFINITE SHUFFLED 
LATTICES IN A STATIC UNIVERSE 

We present in this section results of a numerical study 
of the static toy model above, starting from the simple 
"shuffled lattice" initial conditions. Our aim here in this 
brief presentation is simply to illustrate the qualitative 
similarity of the evolution to that observed in the exactly 
analogous 3 — d model, which has been presented and 
studied in ||2^ (see also 33 ). In particular we wish 



to illustrate that the toy model manifests the "hierarchi- 
cal" and "self-similar" behaviours of the latter, which are 
features also of expanding universe simulations in 3 — d. 
The studies cited above of this case |l|, m, |l§ start 
from different initial conditions and focus on evolution 
on much longer times scales in which the finiteness of the 
box, or initial perturbation, is explicitly important. The 
study in |l8[ [l9| considers, in particular, the regime of 
thermalisation, which is defined only for finite systems. 
On the other hand the dynamics we analyse is very simi- 
lar qualitatively to that studied in the expanding models 
pO| , pT| HJ, [l^ . Indeed the methods of analysis we use 
below are, as in |pT| , the standard ones used in cosmo- 
logical simulations. 



A. Numerical simulations 

In a 3 — d cosmological A^-body simulation, as we have 
noted, the underlying infinite physical system is treated 
numerically using the "replica method", i.e., an infinite, 
but periodic, system is used. The numerical integration 
involves calculating the force on each of the N particles 



by summing over this infinite system. Physical results 
should, of course, not depend on the size of the periodic 
box L. The underlying reason why this is true is that the 
forces on particles converge well in the infinite volume 
limit. More specifically the force on a given particle is 
that due to particles in a finite region about it. The size 
of this latter region is initially of order the mean inter- 
particle separation, but increases as clustering develops 
in the system. As long as the characteristic scale for this 
clustering is small compared to the chosen box size L, 
results are, to a good approximation, independent of L. 
A finite simulation can thus represent well the infinite 
system for a finite amount of time. 

In the 1 — d case we have seen that for a class of per- 
turbed lattices — which are the configurations used as 
initial conditions in cosmological simulations — the force 
is given exactly as a trivial function only of the particle 
displacement. Thus, to simulate numerically the evolu- 
tion of this infinite system, the step in which the force 
is calculated is trivial (rather than involving the approx- 
imation of an infinite sum). The only question which 
arises is how to treat the boundary conditions of the fi- 
nite sub-system of this infinite system which one can sim- 
ulate. Periodic boundary conditions, i.e., particles which 
leave the finite interval on one side enter at the other 
side, are the evident simple choice, as they have advan- 
tage of maintaining (discrete) translational invariance. 
We could, however, easily use other boundary conditions 
(e.g. simply neglecting mass loss, or injecting mass in a 
stochastic manner to compensate average loss), and our 
results should not depend on this choice, just as they 
should not depend on the size of the interval. 

We consider here "shuffled lattice" initial conditions 
and specifically with the PDF for the (independent) dis- 
placements applied to particles from their initial lattice 
sites: 



p{u) 



^ A/2,- 
otherwise . 



-A/2] 



We set the initial velocities to zero. There are in this case 
thus two parameters in the model: the lattice spacing 
I and the amplitude A of the "shuffling" . As we are 
treating the infinite system limit, and the gravitational 
force provides itself no characteristic length scale, there 
is in fact only one relevant parameter characterizing this 
class of initial conditions, which can be taken to be the 
dimensionless ratio y (jus t as in the analogous initial 
conditions in 3 — d, see [P8|). 

Numerically we have simply evolved the particle posi- 
tions as given by Eqs. (|5|) and (|56| ) between crossings. 
The subsequent crossing is determined at each time, and 
the positions and velocities of the crossing particles are 
updated accordingly. For numerical efficiency we have 
implemented the optimized algorithm, using a heap, de- 
scribed in detail in 
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B. Evolution of clustering: visual inspection 

Shown in Fig. ^ are snapshots of the initial conditions 
and evolved configurations a,tt — 2, 5, 8, lOrdyn, for a sys- 
tem with 5000 particles. The plots in the left hand panels 
show the number of particles N{i) in each lattice cell at 
each time, which is proportional to the mass density in 
each cell. Defining the number density contrast as 



Six) 



n{x) — uq 



(74) 



where n{x) = J2iLi Sd{x — Xi) is the microscopic num- 
ber density, the plots represent the evolution of S{x) + 1, 
where the bar indicates an average over the unit lattice 
cell. In the phase space plots, in the right hand panels, 
each point represents simply one particle. 

One sees clearly that the evolution appears to proceed 
in a "bottom-up" manner by the hierachical formation of 
clusters which increase in size in time, starting from the 
smallest scale (of order the lattice spacing). The sense 
in which the system is representative of the evolution of 
an infinite system is manifest in the fact that the system 
does not appear to have a preferred centre — clusters 
form in apparently random locations without sensitivity 
to the boundaries. Indeed we do not follow the evolution 
for longer times than those shown precisely because the 
finite size of the system then becomes important. We 
note that the dynamical features displayed in this regime 
are similar to those seen in the first studies of the 1 — d 
expanding model of M. 



C. Evolution of the Power Spectrum 

In cosmology the primary diagnostic used to charac- 
terize the evolution of clustering in infinite space is the 
power spectrum (PS) (or structure factor) of the particle 
system. Since we consider distributions which are peri- 
odic in an interval of size L, we can write the density 
contrast as a Fourier series 



6{x) — — ex-p{ikx)S{k) , 



(75) 



with k E {{2tt/ L)n,n E Z}. The coefficients S{k) are 
given by 



Sik) 



dx S{x) ex.p{ikx) 



(76) 



The evolution of the PS estimated using an average 
over 500 realisations of our system (with 5000 particles) 
is shown in Fig. |^. 

The evolution is qualitatively like that observed in 3— c? 
simulations (both static and expanding), and the 1 — d 
expanding simulations of ||ll[. At small wavenumbers 
the evolution of the PS shows a simple temporal am- 
plification. This is the behaviour expected from a lin- 
earized treatment (see, e.g., [^) of the equations for a 
self-gravitating fluid, which gives independent evolution 
of each mode of the density field governed precisely by 
Eq.(p3) in the expanding case, and Eq. (54) for the static 
case . 

For the latter case, and vanishing initial velocities, this 
gives 



P{k,t) = P{k,0) COSh\t/Tdyn) ■ 



(78) 



We will see below that this is indeed a very good de- 
scription of the small k evolution. Note that the initial 
PS at small fc is a simple power law P{k) oc k'^, which 
is that of the initial shuffled lattice. The exact form of 
this PS can in fact be derived analytically (see for 
the exact expression, and 35 for a derivation), and 



for the case we are considering takes the form of a sim- 
ple interpolation between this small k behaviour and the 
asymptotic behaviour P{k) = l/ng (characteristic of any 
point distribution). 

The regime in which linear amplification is valid de- 
creases with time, i.e., linear amplification is observed 
in a range k < kNL{t), where k^Lit) is a wave-number 
which decreases as a function of time. This is precisely 
the qualitative behavior one would anticipate as linear 
theory is expected to hold only above a scale which, in 
real space, because of clustering, increases with time. 
At all times, the PS converges at large wave-numbers 
{k > /cjv, where kjsi is the Nyquist frequency) to the 
asymptotic value I/uq. This is simply a reflection of 
the necessary presence of shot noise fluctuations at small 
scales due to the particle nature of the distribution. In 
the intermediate range of fc, i.e. kMhit) < k < fc^r, the 
evolution is quite different than that given by linear the- 
ory. This is the regime of non-linear clustering in which 
the density fluctuations arc large in amplitude. 

One of the important properties of cosmological sim- 
ulations is that, starting from initial conditions with 
PS which are simple power law (for k smaller than the 
Nyquist frequency), the evolution of clustering is, at suffi- 
ciently long times and large scales, "self-similar" (see 
for a discussion and references). By this it is meant that 
the evolution of clustering, above a given spatial scale. 



The PS is then defined as (see e.g. |22 ) 



P(fc) = -(|5(fc)P) 



(77) 



where represents the average over an ensemble of 
realizations of the system. 



Note that, for the 1 — d model the derivation of this result for the 
density contrast follows trivially from the continuity equation. 
The latter gives 5{x) = —du(x)/dx in the limit of small density 
fluctuations <5(x), so that amplification as a function of time of 
S{x) in this limit is simply that of the displacements u(x). 
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FIG. 4: Left panels: the number of points in each lattice cell in the initial conditions (first panel) and at times t = 2,5, 8, lOrdyn', 
Right panels: distribution of particles in single particle phase space at the same times. 
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FIG. 5: Evolution of the PS, averaged over 500 realizations 
of a system with 5000 particles. We have set 1/ = 1. 



is equivalent to an appropriate time-dependent rescaling 
of the length scales. For the PS this relation is written 
conveniently (in 1 — d) as 



kP(k, t) = kRs{t) X P{kRs,ts) , 



(79) 



where Rs{t) is the time dependent rescaling of length^'^, 
normalized by at some arbitrary time tg . Physically this 
behaviour (observed in 3 — ci for a range of initial PS) is 
interpreted as due to the fact that the non-linear cluster- 
ing at a given scale at a sufficiently long time is sensitive 
only to the initial perturbations at larger scales and their 
evolution (and insensitive, notably, to the initial interpar- 
ticle distance, which provides another potential charac- 
teristic length scale in the initial conditions) . Indeed this 
latter evolution (i.e. the regime of linear amplification) 
is, for a power-law initial PS, itself "self-similar" in the 
sense of Eq. (|7^), and the observation of self-similarity 
is simply that this relation extends into the non-linear 
regime. The small k behaviour of the PS (proportional 
to fc^) taken together with the fact that it is amplified 
at small k as given by linear theory then imply that the 
self-similar scaling will be characterised by the function 



Rs{t) = 



COSh(t/ Tdyn 



2/3 



'2t 
exp I -- 



t. 



(80) 



^COSh{ts/Tdyn) J ""'^ \5 Tdyr, 

It is the latter exponential form, for asymptotically large 



Note that this kind of scaling behavior of the PS is found also 
in other statistical physics problems characterized by dynami- 
cal structure formation as for instance coarsening and spinodal 
decomposition in first order phase transitions [p6| . 
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times, which is the relevant one for the self-similar be- 
haviour, as in this limit the reference time is arbitrary. 




R,(t) 



FIG. 6: Evolution of kxP(k. t) as a function of kxRs{t) where 
Rs{t) is as given in Eq. (gO|). Exact self-similar evolution 
would corresponds to the superposition of the curves. 

To assess the validity of this approximation in our sys- 
tem, we show in Fig. |6|the temporal evolution of /c x P{k) 
as a function of the dimensionless parameter k x Rs(t), 
and taking ts = 0. At small k, we see that right from the 
initial time the self-similarity is indeed followed (as the 
rescaled curves are always superimposed at these scales). 
This is simply a check on the result validity of linear 
theory in this regime, as anticipated above. As time pro- 
gresses we see the range of k in which the curves are su- 
perimposed increases, extending further with time into 
the non-linear regime. This is precisely what is observed 
in the analogous 3 — d simulations: as non-linearity devel- 
ops it is characterized by this self-similarity. Note that 
the behavior at asymptotically large k is constrained to 
be proportional to k/no at all times, corresponding to the 
shot noise present in all particle distributions with aver- 
age density uq and which, by definition, does not evolve 
in time (and therefore cannot manifest self-similarity). 



D. Evolution of the Mass Variance 

It is instructive also to characterise the evolution 
through simple real space statistics. In order to distin- 
guish the "non-linear" regime of large fluctuations from 
the "linear" regime of small fluctuations (in which the 
linear fluid theory is expected to be valid), it is useful to 
consider the normalized variance of particle number (or 
mass) in intervals, defined as 



(NHx)) - {N{x)y 



(81) 
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where N{x) is the number of particles in an interval of 
length 2x. 



FIG. 7: Normalized mass variance in intervals of width 2x for 
the different times indicated. As in the previous figures the 
results are for an average over 500 realizations of a system 
with 5000 particles. 



We show in Fig. ^ the temporal evolution of cr^(a;) es- 
timated using an ensemble average over 500 realizations 
of our system (and using periodic boundary conditions 
in the estimation, as in the simulations). As in the case 
of the PS above there are three regimes. At large scales 
we see a simple amplification of the initial functional be- 
haviour, which in this case corresponds to cr'^{x) oc x~'^. 
Note that this behaviour simply corresponds to unnor- 
malized mass fiuctuations independent of scale, which is 
the most rapid decay (proportional to the surface) possi- 
ble in any spatially homogeneous point distribution. At 
small scales, on the other hand, we observe cr^(a;) (x x~^ 
which is the shot noise behaviour intrinsic to any such 
distribution at small scales. The range of scales between 
these two limiting behaviours is that of the non-linear 
clustering resolved non-trivially at this particle density. 
We see that the cross-over to this regime from the linear 
regime occurs approximately where the amplitude of the 
fluctuations is or order unity. 

This behaviour of the variance illustrates very clearly 
the "hierarchical" nature of the clustering, which is 
generic in the evolution of 3 — d simulations starting 
from a very broad class of initial conditions: the initial 
small fiuctuations at a given scale are amplified (accord- 
ing to linear theory) until the fiuctuations in overdense 
regions collapse forming structures. Theoretically such 
behaviour is expected |21 1 for any initial fiuctuations with 
PS with behaviours Pjk) ^ fc" where n < 4, while for 
n > 4 it is expected that the effects of clustering at small 
scales will dominate over that of the linear evolution of 
the very uniform distribution at larger scales. 

It is instructive to probe also in real space the self- 
similar behaviour described using the PS. To do so con- 
sider, following the temporal evolution of scale 
A(q!, t) defined by the relation 



a , 



(82) 
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FIG. 8: Evolution, as a function of time, of the scale at which 
the normalized mass variance has the indicated amplitude. 
The lines correspond to the exactly self-similar behaviour (see 
text). 

where a is a chosen constant. Self-similarity, at a given 
amplitude of the variance, then corresponds to A(q;, t) oc 
Rsit). In Fig. ||we show \{a,t) for the indicated values 
of a, as well as the curves proportional to Rs(t) cor- 
responding to the self-similar behaviour. Such a graph 
illustrates even more clearly than the evolution of the PS 
how self-similarity propagates progessively to the regime 
of larger amplitudes as time goes on, describing in this 
way the clustering further into the non-linear regime. 

While these behaviours are all qualitatively similar to 
those in the 3 — c? case, there are also notable differences. 
For example the maximal amplitude reached by the vari- 
ance in the range of times explored is much smaller than 
that in 3 — d — of order only a few compared to a hun- 
dred or more in the latter case|^. This is probably 
indicative of a much weaker non-linear clustering in the 
former case, associated with the smoother behaviour of 
the 1 — d force at small scales. This difference is also 
reflected in the behaviour of the PS, which flattens much 
more rapidly to its asymptotic Poisson value in the 1 — d 
case — at around the Nyquist frequency rather than at 
considerably larger wavenumbers in the 3 — d case. 



VI. SUMMARY AND DISCUSSION 

We have revisited in this article a basic question con- 
cerning the definition of the gravitational force in 1 — d 
infinite point distributions. Previous definitions of this 
quantity in the literature have required the assump- 
tion of the existence of a special point (centre) in the 
distribution, i.e., explicit global breaking of statistical 
translational invariance which is typically a feature of 
the infinite distributions one instead wishes to study. 
We have noted that the problem, associated with the 
non-converging surface fluctuations in such distributions. 



may be solved by employing a definition using a smooth 
screening which is sent to zero at the end of the calcu- 
lation. We have then shown explicitly that this leads to 
a well defined force for a specific class of infinite per- 
turbed lattices — those subject to perturbations of fi- 
nite variance which do not make particles cross. In this 
case, when the mean displacement of particles is also as- 
sumed to vanish, the force on each particle take a unique 
value which is simply proportional to its own displace- 
ment from its lattice site. We note that we have assumed 
also that variance of the displacement fields is finite, 
which restricts to initial density fiuctuations which have 
a sufficiently rapidly decaying power spectrum at small 
wavenumber (specifically, such that P{k — s-) fc" where 
n > 1, analogous to the same condition with n > — 1 in 
3-d). 

We have then discussed different dynamical toy mod- 
els which incorporate this definition of the force — the 
simple conservative Newtonian dynamics and then one 
which incorporates a damping term mimicking the effect 
of 3— d expansion. Since the crossing of particles is equiv- 
alent, up to labels, to elastic collisions with exchange of 
velocities, the configurations generated by such dynam- 
ics, at any finite time, are always in the class of infinite 
perturbed lattices for which the force is defined (provided 
such a configuration is the initial condition). This is the 
case because, at any finite time, collisions/crossings may 
only correlate particles up to a finite distance, and the 
correlation properties of displacements at asymptotically 
large separations therefore always obey the required con- 
ditions. The equations of motion are then simply those 
of an infinite set of inverted harmonic oscillators (with 
damping in the expanding case) with centres on the orig- 
inal lattice sites, and which bounce elastically when they 
collide. In this context we have also discussed in detail 
the different formulations of these models in the previous 
literature. We have then presented some results for the 
development of clustering in the simplest model, for the 
simplest class of initial conditions. We have underlined 
the similarity of the results to those which have been ob- 
tained for the analogous system in 3 — d, which is itself 
a simplified model for the full cosmological model. The 
physical meaning of the prescription adopted for the force 
is also very clearly illustrated by these simulations: the 
dynamics we observe in this infinite system limit is sim- 
ply that which would occur in a system with a screened 
gravitational force, in the regime in which the size of the 
structures is much less than the scale of the screening. 

A few additional remarks on our prescription for the 
force, and the class of point distributions we have con- 
sidered, are appropriate: 

• While we have emphasized that our definition of the 
force does not require explicit breaking of transla- 
tional invariance, we have in fact shown it to give 
a well defined answer only for a class of point pro- 
cesses which have discrete statistical translational 
invariance (by lattice vectors) rather than full sta- 
tistical translational symmetry. We will consider in 
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forthcoming work the more general question of the 
general conditions on point processes for the force 
to be defined. In this respect we note that it has, in 
fact, been shown in that any point process in 
1 — d with a variance of mass in an interval which 
is bounded necessarily breaks continuous statistical 
translational invariance. While we have not shown 
here that the definedness of the 1 — d force requires, 
in general, such boundedness of the variance, it is 
straightforward to show that this boundedness in- 
deed characterizes the class of perturbed lattices 
for which we have found the force to be defined. 

• We note that the results of |^ in fact generalize 
an earlier result (see and references therein), 
that the thermodynamic equilibria of the 1 — d one 
component plasma (OCP, or"jellium" model) break 
translational invariance. This model is in fact just 
the same system we are considering here up to the 
sign of the force, with the difference that the pres- 
ence of a physical neutralizing negative background 
is specified. We therefore do not have, a priori, the 
freedom to use the regularization of the problem in 
the infinite system limit which we have exploited 
here. We note, however, that if we were to do so, 
i.e., define the model in the thermodynamic limit 
as the zero screening limit of a screened Coulomb 
interaction in 1 — d, we obtain, for the class of per- 
turbed lattices we have determined, a quite trivial 
model in which all particles are simple harmonic 
oscillators about their lattice sites, which bounce 
elastically when they collide. It would be interest- 
ing to investigate further whether such a formula- 
tion of the OCP is in fact equivalent to the usual 



• Even if the force itself is not defined, with our pre- 
scription, for a given point process in the infinite 
volume limit, it may still be possible to construct 
toy models of clustering in an infinite space, with 
only the weaker condition that differences in forces 
at a finite separation be defined. This is related 
also to the restriction on the variance of the dis- 
placements we imposed, which we have noted cor- 
responds to a restriction on the small wavenumber 
behaviour of the power spectra for which the force 
is defined. If one considers differences in the forces 
rather than forces, it is simple to show that any 
power spectrum of density fluctuations which is in- 
tegrable at small wavenumber can be treated. We 
will discuss this issue, and the more general con- 
ditions for definition of the force, in forthcoming 
work. 

• In cosmological simulations initial conditions are 
prepared by applying stochastic perturbations to 
a perfect lattice (see e.g. Is^, |40|). These per- 



turbations usually have a correlated Gaussian joint 
PDF, which means that in the 1 — d analogy the 
no-crossing condition we have required cannot be 
satisfied. In practice, however, the initial ampli- 
tudes are always chosen sufficiently small so that 
no such crossing occurs in the finite simulation box, 
and indeed the approximation in which the algo- 
rithm for the initial conditions is valid corresponds 
to this limit. To establish rigorously that our re- 
sults can be extended to incorporate this case would 
require a generalisation of the calculation we have 
given which incorporates the contribution to the 
force PDF from crossings. This will lead to a non- 
zero variance in the PDF, but one would expect 
that its effect should indeed be negligible in the 
limit that the typical displacement (measured by 
the displacement variance) is small compared to the 
lattice spacing. 

In future work we aim to exploit further this model, 
and the expanding variant, as toy models for the cos- 
mological problem. Despite the previous works in the 
literature which have explored various versions of these 
models in the regime in which the analogy to 3 — d cos- 
mological simulations is most direct, and found various 
simple behaviours in the clustering — notably the stud- 
ies of 0, |l], |l|, — they have led so far to httle of 
the analytical insight one might hope to gain from a toy 
model. For example, although fractal behaviours have 
been documented in numerical simulations by ||l^, |lj, , 
the relevant exponents remain unexplained. The very 
simple formulation of the models we have given may help 
in this respect. 

These toy models may also be useful in understanding 
other aspects of numerical simulations in 3— d, such as the 
question of the importance of discreteness effects: in cos- 
mological simulations the aim is to reproduce as closely as 
possible through a particle simulation the Vlasov-Poisson 
limit, which corresponds to an appropriate infinite par- 
ticle number limit. The corrections which arise due to 



finite particle number are poorly understood (see |41| for 
a detailed discussion and references), leading to the ab- 
sence of control on the associated errors in theoretical 
predictions. Although there are evidently important dif- 
ferences, the problem can be formulated and studied in 
the greatly simplified 1 — d context. This requires firstly 
a clear formulation of the Vlasov-Poisson limit, which so 
far has been given rigorously only for finite systems p2[ . 
This then allows one to define an appropriate numerical 
extrapolation which should be used to study convergence. 
In 1 — d there is the interesting possibility also of perform- 
ing directly simulations of the Vlasov-Poisson system for 
comparison. 

We thank D. Fanelh, M. Kiesshng, B. Marcos, B. Miller 
and F. Sylos Labini for useful discussions and exchanges, 
and B. Jancovici for providing us with references ||37|, p8[. 
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